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^ ' Abstract 
> ■ 

i It is shown that the renormahzation group (RG) method for global analysis can be 

I formulated in the context of the classical theory of envelopes: Several examples from 

' partial differential equations are analyzed. The amplitude equations which are usually 

^n , derived by the reductive perturbation theory are shown to be naturally derived as the 

^ I equations describing the envelopes of the local solutions obtained in the perturbation 

Q ■ theory. 

» 

-I— > 

1 Introduction 

> 

■ Renormalization group (RG) equations|l|] appear in various fields of science. In quantum 
field theory]^], the RG equation improves results obtained in the perturbation theory. In 
recent years, the improvement of the effective potential 0, 0, |[ has acquired a renewed 
interest p. The RG equation has also a remarkable success in statistical physics especially 
in the critical phenomena 0. One may also note that there is another successful theory of 
the critical phenomena called coherent anomaly method (CAM)0; CAM utilizes envelopes 
of susceptibilities and other thermodynamical quantities as a function of temperature. It 
is well known that Feigenbaum applied RG equation to deduce a universality of some 
chaotic phenomena. 

Recently, Illinois group [§, ^ have shown that the RG equation can be also used for non- 
quantum mechanical problems: They proposed to use the RG equation to get an asymptotic 
behavior of solutions of differential equations including ones of singular and reductive per- 
turbation problems in a unified way. Mathematically, the RG equation is used to improve 
the global behavior of the local solutions which were obtained in the perturbation theory. 
This fact suggests that the RG method may be formulated in a purely mathematical way 
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without recourse to the notion of the RG. In the previous paper [|rT[|, we showed that the 



RG method can be formulated in the context of the classical theory of envelopes fl^: We 
pointed out that the RG equation is nothing but the envelope equation, and gave a proof 
why the RG equation can give a globally improved solution to ordinary differential equa- 
tions(ODE's). In fact, if a family of the curves {C^}r in the x-y plane is represented by the 
equation F{x,y,T) = 0, the equation G{x,y) = representing the envelope E is given by 
eliminating the parameter r from the equation 



dF{x,y,T) 



dr 



0. (1.1) 



Here we have chosen the parameter to be the x-coordinate of the point of the tangency 
of a curve and the envelope E. The relevance of the envelope equation Eq.(l.l) and 
the (Gell-Mann-Low type) RG equations is evident; the parameter r corresponds to the 
renormalization point. Thus one would also readily recognize that improving the effective 
potential in quantum field theory is constructing the envelope of the effective potentials with 
the renormalization point varied. |Tl|] 



The purpose of the present paper is to show that the formulation can be naturally ex- 
tended to a class of partial differential equations (PDE's); the PDE's dealt in the present 
work include a dissipative nonlinear hyperbolic equation, one- and two-dimensional Swift- 
Hohenberg equations, damped Kuramoto-Shivashinsky equation[|l^] and Barlenblatt equation [p!^ 



Most of these examples were examined by the Illinois group P, ^, |T0| . Therefore there will be 
some overlaps in the exposition with theirs. It is inevitable, however, because our purpose is 
to give a purely mathematical formulation of their method without recourse to the notion of 
the renormalization group. We will show that the amplitude equations usually obtained by 
the reductive perturbation methods are given naturally as envelope equations. In the course 
of the formulation, relevance of the characteristic manifold as the branch stripmanifold will 
be indicated to constructing the global solutions. In Appendix A, we give a short summary 
of the classical theory of envelope surfaces; for PDE's, envelope surfaces may be more rel- 
evant in some cases than envelope curves discussed in In Appendix B, we examine 
Barenblatt's equation in our approach; for this equation, only the anomalous exponent of 
the long-time behavior of the solution is given. 



2 Dissipative nonlinear hyperbolic equation 



We first consider the following slightly dissipative nonlinear hyperbolic equation [^] to apply 
the envelope theoryfl to construct a global solution: 

+ \{u)^^ — rj^"^^ (2 1) 

dt dx dx'^ ' 



^See Appendix A and for classical theory of envelopes. 
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where X{u) is a sufficiently smooth function of u, and 77 is a positive constant. Following 
we consider a small amplitude wave in the background of constant solution uq] 

u = Uq + eui + e^M2 + • • • . (2.2) 

Here e denotes the amplitude of the wave. \{u) is expanded as 

A(M) = Ao + eA'(Mo)Mi + ---, (2.3) 

where Aq = A(mo). We further assume that the dissipation is weak in the sense ~ e; we 
write as ?7 = fie, where fi = 0(1). Then equating the coefficients of e" (n = 0, 1,2, ...), we 
have 

dtUi + Ao^^Mi = 0, 

dtU2 + XodxU2 = -X'QUidxUi + iidlui, (2.4) 
and so on. Here, Aq = A'(uo). 

It is now convenient to introduce the new variable ^ = x — XqI, the characteristic direction 
of the unperturbed equation; we describe the solutions in terms of (^,t). Then 

dtUi = 0, 

dtU2 = -X^uid^ui + fid'^ui, (2.5) 

and so on. One readily gets 

wi(e,t) = Fo(o, 

U2{U) = {t-to){-X',F,{Od^Fo{0 + f^dlFo{0), (2.6) 
where -fo(0 is an arbitrary function of ^. We note that U2 is a secular term.l 

We have thus an approximate solution to Eq. (2.1) 

n(e, t; to) = uo + eFo(0 + e\t ~ to)(-A[,Fo(e)aeFo(0 + /ia|Fo(0) + 0{e^), (2.7) 

where we have made it explicit that u is dependent on an arbitrary time to- We stress that 
this (approximate) solution is only valid for t around to. 

Now, geometrically speaking, we have a family of surfaces S'fp represented by u{C,,t]to) 
with a parameter to. Let us obtain the envelope surface E of this family of surfaces, following 
the classical theory of envelopes given in Appendix A: we represent E by Mg(^,t). 

We ffist note that -Fo(0 t^^Y be functionally dependent on the arbitrary initial time to: 

M(e, t; to) =uo + eFo(e, to) + e\t - to)(-A[,Fo(e, to)d^Fo{^, to) + /i9|Fo(e, to)) + 0{e^). (2.8) 

^One might have added another arbitrary function Go{£,)oi ^ to U2. However, the effect of Go could be 
renormahzed away to Fq. 
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It is natural to set the tangent curve Ctg of E and Stg to lie along t = to-line, which is parallel 
to the characteristic direction (^-direction): 

Ct,: t = to, u = u{^,to,to), e = e (2.9) 
Then the envelope equation reads 

— = 0, and to = t, (2.10) 

OlQ 

which leads to 

dtF{^, t) + eX'iuo)Fd^F = r/a|F. (2.11) 

This is Burgers equation which is usually obtained in the reductive perturbation theory. 
Then 

«^(e, t) = n(e, t; t) = uo + eF(e, t). (2.12) 



Now one may wonder if Uj^{^,t) satisfies Eq. (2.1), although M(^,t;to) does up to O(e^) 
for any to- We remark that the question is not trivial because M^(^,t) = u{^,t]t). We shall 
show that the answer is yes. In fact, the time derivative of Uj^{^, t) at t = Vto coincides with 
du{^, t; to)/dt at t = to; 



du^{i,t) 



dt 



du{^,t;to) 



t=to 



dt 

du{^,t;to) 



du{^,to;t'Q) 



dt 



t=to 



t=to 



dt'n 



t'o=to 



(2.13) 



where Eq. (2.10) has been used. Furthermore, needless to say, Uj^{C,,to) = u{C,,to,to). Thus, 
at t = Vto, 



^^^+A(«^)?^ = r,?% + 0(a. 



dt 



dx 



dx^ 



(2.14) 



Thus one sees that our envelope function satisfies the Eq.(2.1) uniformly for t in the 
global range. 



A comment is in order here: one may say that to get the approximate but global solution 
Uj^ from the local solution Eq.(2.8), we have utilized the fact that the uniqness of the solution 
of the Cauchy problem is violated when the initial data are given along the characteristic 
manifold (curve); see Appendix A for characteristic manifolds(curves). 



3 One-dimensional Swift-Hohenberg equation 



In this section, we deal with the one-dimensional Swift-Hohenberg equation 

Liu = eu- u\ Li = dt + {dl + ^f, 
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(3.1) 



where e is a small parameter. We shall show that the envelope of local solutions of Eq. (3.1) 
satisfies the time-dependent Ginzburg-Landau equation. 

Following we scale u as 

u = ye0. (3.2) 

Then (f) satisfies 

L,0 = e(0-03). (3.3) 
We solve this equation in the perturbation theory, expanding (p as (f) = (f)Q + e(j)i + we have 

Li0o = 0, = 00 - 00, (3-4) 
and so on. As the 0-th order solution, we take 

0o(a;,t) = Ae^'^^ + c.c, (3.5) 
where c.c. denotes the complex conjugate. Then the Ist-order equation reads 

Li0i = Ae"''' - A=^e=^*^" + c.c, (3.6) 

where 

A = A{l-3\A\''). (3.7) 
A special solution to Eq. (3.6) is found to be0 

01 = A{fi,it - to) - ^(^' - 4)}e''^ + ^,e'^'^ + c.c, (3.8) 

with fii + fi2 = 1. We note that the secular terms have appeared in 0i because of the first 
term in r.h.s. of Eq. (3.6). 

Now one may say that we have a family of surfaces St^xo represented by 

0(x,t;xo,to) = [Ae''^ + eA{fi,{t-to)--^ix'-xl)}e''^ + e-^e'^''^] 

+C.C., (3.9) 

parametrized by to and Xq. Let us obtain the envelope of the family of the surfaces in two 
steps, by assuming that the amplitude A is dependent on xq and to. First, fixing xo, we 
obtain the envelope Ei of the surfaces with to being the parameter. The resultant envelope 
has the parameter xo] then we obtain the envelope of the family of the surfaces Ei. 

The first step is achieved by setting 

^ = 0, with to=^- (3.10) 

dto 



This is the condition to get the envelope that has the common tangent curve with Sx^tg along 
a;-direction. Thus we have the equation for A{xo,t); 



dtA = fi,eA + 0{e^) 
Then the envelope Ei is represented by 

(l)Ei{x,t;xo) = (p{x,t;xo,to = t), 



{3.11] 



\A{x^,t)e^'^ - eA{x^,t)r±-(x' - x^Y'^ + e 



+ C.C.. 



(3.12) 



The envelope E of the family of the surfaces given by Eq. (3.12) is obtained as follows; 



dxv, 



0, with xo = X, 



which leads to 



X 



dxA[x,t) = -^^eA^ + Oie^) 



(3.13) 



(3.14) 



Here we have utilized the fact that d^A is 0(e). Differentiating this equation with respect 
to X, we get 



dlA{x,t) = -e^^A + 0{e'). 

Here we note again that d^A is 0(e). 

Now combining Eq.'s (3.10) and (3.14), we have 

dtA = 4k^dlA + eA{l - 3\A\^), 
up to O(e^). The envelope is represented by 

0ij(x,t) = (f)Ei{x,t;xo = x), 

A{x,t)e^'^ + e^^;^e'''^ 



+ c.c. 



(3.15) 



(3.16) 



(3.17) 



A comment is in order here: To derive Eq.(3.15), we started from the first-order differen- 
tial equation Eq.(3.13) with respect to Xq. One could start from the second-order differential 
equation as is done in 



dXn 



0, 



(3.18) 



X(}=X 



which gives Eq.(3.15) as an exact relation without the remainder of O(e^). We remark that 
the geometrical meaning of Eq.(3.18) is not clear in contrast to Eq.(3.13). In our formulation, 
Eq.(3.18) is derived as an approximate relation, as is shown below. 
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Now one may ask the question as to the relation of (f)E{x,t) given in Eq. (3.17) and 
the original equation Eq.(3.1): The answer is that (j)E{x,t) satisfies Eq.(3.1) up to O(e^) 
uniformly in the global domain due to the very envelope conditions Eq.(3.10) and Eq.(3.13). 
In fact, for \/t = to, 



dt 



t=to 



dt 



t=to 



dto 



t=to 



dt 



(3.19) 



t=to 



on account of Eq.(3.10). Similary, for Vx = xq, one can easily verify that 



92 



Q2 



X=X(} 



dx"^ 



(3.20) 



x=xo 



on account of Eq.'s (3.12), (3.13) and (3.18). For instance, d'^(j)E/dxdxo\x=xo = O(e^). 



4 Damped Kuramoto-Shivashinsky equation 



In this section, we deal with the (one-dimensional) damped Kuramoto-Shivashinsky equa- 
tion, which is given by 

Liu = eu — udxU, (4.1) 
where Li is defined in Eq.(3.1). This equation is not examined by the Illinois group. 

By scaling u as 

u = V~e<P, (4.2) 



we 



havei 



Lie 



.1/2, 



+ e(h. 



(4.3) 



We first try to solve this equation by the perturbation theory with the expansion 

(j) = (j)o + e^/^(j)i + e(f)i + .... 
The equations for {/)o, <^2--- are found to be 

Ll0O = 0, Ll0l = -0o(?x0O, -^^102 = 00 - (0O<9x01 + 0lC^x0o), 

and so on. 



(4.4) 



(4.5) 



As the 0-th order solution, we take 



Ae*^"" + c.c. 



(4.6) 



^If one scales u as u — ecf), one will get a different equation from that given below. 



where A is a constant. Then 0i is found to be 



01 = - 

Thus the equation for 02 reads 

Li02 = A{1 - 



\A\'^ \Jkx 



(4.7) 



(4.8) 



which is in a similar form with Eq.(3.6). One easily gets for a special solution to this equation 

re™ + C.C.. (4.9) 



Thus we reach the solution in the perturbation theory up to 0(e^/^), 



(j){x,t;xo,to) 



Ae 



A^ 
192F' 



^2 / 2 ^2\'[ikx 



8k^ 



Jiikx 



+ C.C.. 



(4.10) 



Now we have a family of surfaces Sx^ta in x-t-(j) plane represented by (f){x,t; Xo,to) with 
Xo and to being the parameters. We repeat the procedure of the previous section to obtain 
the envelope E of the family of the surfaces: We first note that A may depend on xq and 
to, i.e., A — A{xo, to). Then fixing xq, we first obtain the envelope of Sxoto with to being the 
parameters. The envelope Ei is obtained by 



0, with to — t, 



(4.11) 



where we have assumed that the tangent curve (the characteristic curve) is along x-direction 
by setting to = t. The above equation leads to 



dtA{xo,t) = ef,,A{l - \^) + 0{e'/'). 

With A satisfying this equation, the envelope Ei is represented by 

4>Ei{x,t]Xo) = (j){xo,to;xo,to = t) 

Ae'^- - te^^^e^^^- - eA{l - ^)-^(x2 - xl)e'^- 
9k^ ^ 9A;2^8A;2^ °^ 



(4.12) 



— e 



A^ 
192F' 



+ C.C. 



(4.13) 



This function can be regarded as representing a family of surfaces with xo being the param- 
eter. The envelope E of this family of surfaces are obtained by setting 



dxo 



0, with Xq = X, 



(4.14) 
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which leads to 



(4.15) 



Here we have utihzed the fact that d^A ~ 0(e). Further differentiating with respect to 
one has 



Combining Eq.'s (4.12) and (4.16), one sees that the amphtude satisfies equation 

{d,-Ak'dl)A^eA{l-\^), 

up to O(e^). 

With this amphtude, the envelope E is given by 

(j)E{x,t) = (j){x,t]Xo = X), 

192A;6 



(4.16) 



(4.17) 



(4.18) 



A couple of comments are in order here: 

(1) Eq.(4.16) could be obtained as an exact relation by imposing the condition that 



dxi 



0. 



(4.19) 



xo=x 



We must note, however, that the geometrical meaning of this condition is not clear. In 
contrast, in our formulation, the second derivative is evaluated to be 



9^0 



dxl 



(4.20) 



Xq=X 



(2) As was done in the preceding section, one can easily show that (f)E{x,t) satisfies the 
original equation Eq.(4.1) up to O(e^) but uniformly in the global domain. 



5 Two-dimensional Swift-Hohenberg equation 



In this section, we deal with the two-dimensional Swift-Hohenberg equation given by 

L2U ^eu- 4 = 5* + {dl + + k^f. (5.1) 
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Scaling u as 

u = A/e0, 

we have 

We first solve this equation in the naive perturbation theory expanding = 0o+e0i+e^02+---; 

1/2(^0 = 0, L20„ = - (n = l,2, ...). 

If we assume the roll solution along y-axis for the zero-th order equation, 

00 = Ae^'^^ + C.C., 
then the solution up to O(e^) is found to be [l^ 



(5.2) 



(5.3) 



(5.4) 



(5.5) 



A + e{/ii(t - to) - fJ'2 



+ At3 



8ik 



+ fJ'4. 



y - Vo 

4! 



}A 



Akx 



+C.C., with A = A{l~^A\^), 



(5.6) 



where I]i=i~4/^t = 1 and Xo,l/o and to are arbitrary constants. Here we have omitted the 
terms that do not give rise to secular terms. 

Now one may regard that we have a family of "surfaces" in the x-y-t-cj) space with xo, yo 
and to being the parameters. Let us obtain the envelope of this "surfaces" in the three steps 
by noting that A may be dependent on Xo,|/o and to: First we regard that only to is the 
parameter of the "surfaces" with both xq and yo fixed. The condition reads 



atn 



0, tn = t. 



(5.7) 



This condition may be also regarded as the one for the envelope curve of a family of curves 
in t-0 plane with x and y being fixed. The condition leads to 



dtA{xo,yo,t) = efiiA + 0{e^). 
Inserting this solution, we have the envelope 

4>Ei{x,y,t;xo,yo) = (f){x,y,t; xo,yo,to = t) 



{5.t 



A{xo,yo,t) + e{-^2 

+ C.C., 



8P 



xy'^ - x^yl 
8ik 



4 4 

y -Vo 

4! 



}A 



Akx 



(5.9) 



which we may regard as a family of curves in x-(j) plane with xq being the parameter where 
y and yo are fixed. The envelope E2 for this family of curves is given by setting 



dXn 



0, with xo = X, 



(5.10) 
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which leads to 



^4fc2 8ik 
Accordingly, the envelope is represented by 

(l)E^{x,y,t;yo) = 4)E^{x,y,t; Xq = x,yo) 

Akx 



(5.ii: 



A(x,2/o,t)e^'^^ + e{/i3X^^— ^ + /i4^^7T^MeH +C.C., (5.12) 



%ik ^ 4! 

which is regarded as representing a family of curves in y-cj) plane. The envelope of is 
obtained as usual by setting 

d(pE2 



dyo 



0, with yo = y, 



which leads to 



dyA{x, y, t) = e{-^xy + ^^y')A + 0{e'). 



(5.13) 
(5.14) 



Differentiating Eq.'s (5.11) and (5.14), we have 

dlA{x,y,t) : 



d^dlA{x,y,t) - 
d^A{x,y,t) = eiXiA + 0{e 



-e-^A + 0{e') 



where we have utilized the fact that 

d^A ~ 0(e), dyA ~ 0(e). 
Combining these equation together with Eq.(5.8), we finally reach 
dtA - Ak'^dlA + Aikd^dlA + d^A = eA{l - 3\A\^) 

up to 0(e2). 

With this amplitude, the envelope is given by 

(f)E{x,y,t) = (f)E,{x,y,t;yo = y), 
= A{x,y,t)e^'''' + c.c. 



(5.15) 

(5.16) 
(5.17) 

(5.18) 

(5.19) 



(5.20) 



A few comments are in order here: 



(1) By imposing that the any order of the differentiation of with respect to Xq and yo 
should vanish, one could get Eq.'s (5.15-17) as exact relations although the geometrical 
meaning of these conditions are unclear. In our formulation, instead, the following are 
derived as approximate relations 



d'^Xo 



0(e^ 



E2 



xo=x 



dxd^yo 



0(e' 



■>E2 



yo=y 



d^yo 



O(e^). 



(5.21) 



yo=y 



(2) It can be shown that 0^; satisfies Eq.(5.1) up to O(e^) but uniformly in the global 
domain owing to the above relations together with Eq. (5.8). 
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6 A brief summary and concluding remarks 



We have shown in the present paper that the RG method for global analysis can be formu- 
lated for partial differential equations, too. An interesting equation which is treated by the 
Illinois groupP but left untouched in the text is Barenblatt's equation. A complete global 
solution to this equation is not given in the RG method but only the anomalous exponent 
for the long-time behavior is obtained. For completeness, we show that the same result 
for the anomalous exponent can be obatined in our envelope theory in Appendix B. Thus 
together with our previous paper [|ll|] where ordinary differential equations are discussed, we 
have shown that almost all the classes of problems treated in the RG method by the Illinois 



group P, P, 10| are nicely formulated on the basis of the classical theory of envelopes without 



recourse to the notion of the renormalization group. Actually, this is natural because the 
resulting equations indeed describe the amplitudes of the nonlinear waves as given by the 
solutions of the nonlinear equations. 



It is already indicated P, ^ that the renormalizability is equivalent to the solvability of 
equationsJTO . It would be interesting to apply the theory developed here to systems of 
equationsQ and see possible geometrical meaning of the solvability condition. 



In deriving the whole envelopes of local solutions, we have taken a multi-step approach 
in sections 3, 4 and 5. It is interesting that a kind of multi-step approach is also proposed 
for improving the effective potentials in quantum field theories with multi-scales; see the 
paper by Fordp. It may imply that our approach that identifies the RG equation with the 
envelope equation naturally leads to the effective potentials as given by Ford when applied 
to quantum field theories. |T8|] 
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Appendix A 



In this appendix, we give a short review of the classical theory of envelopes 1 12]. We 
shall first consider envelopes in three-dimensional space, i.e., envelope surfaces. Then the 
extension to higher dimensional cases will be briefly described.^ 

Let {Sr}T be a family of surfaces parametrized by r in the x-y-u space; here Sr is 
represented by the equation 

^x,y,u-T) = 0. (A.l) 
We suppose that {Sr}T has the envelope E, which is represented by the equation 

^{x,y,u) = 0. (A.2) 
The problem is to obtain y, u) from <l>(x, u; r). 

Let E and Sj- has a common tangent plane on a common curve C^; with r being varied, 
{Cr}T forms E. Cr is called a characteristic curve. The necessary condition for 5*,- to have an 
envelope E is given as follows. Let Cr be represented by x = x(cr, t), y = y{cr, t), u = u{cr, r), 
with a being a parameter. We assume that 

rank f^'^^'^'''^') =2, (A.3) 

\Xr yr Ur J 

in order for {C^},- to give a non-degenerate surface E. Now, since Cr is on Sr, 

T),y{cr, T),u{a, t);t) = 0. Differentiating this equation with respect to r, one has 

^^Xr + ^yyr + ^uUr + $r = 0, (A.4) 

where = d^/dx and so on. On the other hand, since {xr, yr, Ut) = is a tangent vector 
both of E and Sr, and {^x, ^y, ^u) = n,- is a normal vector of Sr, 

■ = ^xXr + %yr + KUr = 0. (A.5) 

Combining Eq.'s (A.4) and (A.5), one has 

^ d<^{x,y,u-r) ^ 

OT 

as a necessary condition for {Sr} to have an envelope E. Thus, solving Eq. (A. 6) for 
r = t{x, y, u), one finally gets 

^(x, y, u) = $(s, y, u; t{x, y, u)) = 0. (A.7) 



Accordingly, the characteristic curve is given by the conditions; 

\l/(x, y, m) = and r(x, -u) = const. (A. 8) 



In 1 11 , a review is given on how to construct envelope curves. 
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Note that the second equation gives a constraint on x, y and u. 

We remark that the sufficient conditions for the existence of an envelope is supplemented 
by [0 



$rr = 0. (A.9) 

When Sr is given by m = f{x,y; r), the condition Eq. (A. 6) is reduced to 

M^ = 0, (A.10) 

or 

which gives r as a function of x and y. Thus we get for the envelope 

The characteristic curve is accordingly given by 

u = ^^{x,y) (A.12) 

T(x,3/)=const. 

As an example, let 

u = ip{x, y, t) = exp(r?/)(l + y{x — r)) + exp(— x). (A. 13) 

Then the equation dip /Or = gives 

T = x, (A.14) 
except at y 7^ 0. Thus the envelope is given by 

u = exp{xy) + exp(— x). (A. 15) 

The characteristic curve C,- is given by 

X = r = constant = xq, and u = exp(xo y) + exp(— xq). (A. 16) 

We remark that d'^ip/dr'^ ^ provided that y ^ 0. 

The theory can be extended to the envelope of a family of hyper-surfaces {Sr}T in the {n+ 
l)-dimensional space R"-"*"^. Let {Sr}T be represented by the equation ^{xi,X2, ■ ■ ■ , Xn, u; r) = 
with a parameter set r = (ri, r2, ■ ■ ■ , t„_i). Then the function representing the envelope is 
given by eliminating parameters Xj (i = 1, 2, ■ ■ ■ , n — 1) from 

(9$ 

$(xi,X2,---,a;„,M;r) = 0, 7^ = 0, = 1, 2, ■ ■ ■ , n - 1). (A.17) 

OTi 

A few comments are in order here: 
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(i) The envelope of a family of surfaces has usually an improved global nature compared 
with the surfaces in the family. So it is natural that the theory of envelopes may have some 
power for global analysis. 

(ii) It should be stressed here that Eq.(A.6), which we call the envelope equation, has the 
same form as renormalization group (RG) equations. Rather, it should be said conversely; 



RG equations in general are envelope equations []TT|. This is the reason why RG equations 
"improve" things; one should note that the improvement by an RG equation usually means 
that a function with a better global nature is constructed from functions with a local nature 
only valid around the renormalization point. 

(iii) Let u = (p{x, y; a, r) be a complete solution of a PDE 

F{x,y,u,p,q) = 0, {p = du/dx, q = du/dy), (A. 18) 

where a and r are constant. Assuming a functional dependence of a on r, i.e., a = A{t), 
one can obtain the envelope as follows; 

du dA du du ^ ,^ 
dr dr da dr 

from which r is given as a function of x and y, namely, r = T{x,y). Then one sees that 
the envelope function u = (pj^{x,y; [A]) = {p{x,y, A{T{x,y)),T{x,y)) is the general solution 
to F{x,y,u,p,q) = 0. 

Appendix B 

In this Appendix, we examine a non-linear partial differential equation which has an 
anomalous exponent in the long-time behavior of its solutions; 

(d,-^dl)u{t,x)=9{-dtu)dlu. (B.l) 

This is called Barlenblatt's equation. |ll3|, H This equation is a nonlinear PDE with the step 
function in r.h.s. Without this term (e = 0), the equation is a simple diffusion equation and 
has a self-similar form after a long time; u{x, t) ~ t-^l^f{x^/t). However, with e ^ 0, the 
long-time behavior is found to be 

2 

(B.2) 

The anomalous exponent a can be interpreted as an anomalous dimension P]. The appear- 
ance of the anomalous dimension is due to the fact that the scale / with a spatial dimension 
which characterizes the initial distribution of u{x, 0) can not be neglected even in the long- 
time limit t — > cxD in contrast to the case of the simple diffusion equation JTBI. We are 



interested in determining the anomalous dimension a in the perturbation theory. 

To solve the problem, we shall follow for a while. First we convert the equation to an 
integral equation; 



u{x,t) = / dyG{x - y,t)u{y,0) 

J — oo 

g rt roc 

+^ ds dyG{x -y,t- s)e[-dsuiy, s)]d^u, (B.3) 

/ Jo J-oo " 
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where 

G(x,t) = ^^e-^ (B.4) 

is the Green's function! We expand u as 

u{x, t) = uo{x, t) + eui(x, t) + ■ ■ ■ . (B.5) 
As an initial condition, we take 

m(x,0) = Mo(a;,0) = -^==e"5i^, (B.6) 

where tuq and / are parameters. Then one finds that 

t) = ^Le-l^{l - In U + 0(e2) + 0(/Vt, e), (B.7) 



where we have retained in ui only the term which may contribute to the anomalous dimen- 
sion. 

We renormalize mo at t = to as follows; 

m(to) = ^(to, /)mo, Z{to, = 1- In ^, (B.8) 

accordingly, 



M(0,to) =m(to)/V2^, (B.9) 

and then 

«(a:, t- to) = ^e-^{l - -L= In f } + ©(e^) + 0(1^1, e). (B.IO) 

Here we have made the to dependence of u explicit. 

We now leave the line of argument of Goldenfeld et al|]^. We have a family of surfaces 
represented by u{x,t;to); the surfaces are parametrized by to- Let us obtain the envelope 
E of this family of surfaces. We can suppose that the tangent curve Ct^ (the characteristic 
curve) is along t = to: 

t = to, u = u{x, to; to). (B.ll) 



^If the Green's function is defined as usual by 

1 



where 9{t) is the step function. 
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Now according to the theory given in Appendix A, the envelope is given by the following 
condition; 

Oil 

— = with to = t. (B.12) 

This equation is reduced to an equation for m, 

m I 

m + a — = 0, a = e/v27re, (B.13) 

hence 

m{t) = mr" (B.14) 

with rh being a constant number. Thus the envelope is given by 

u^{x, t) = u{x, t- 1) = mt-^^/^+'''^F{x^/t), (B.15) 

hence the anomalous exponent reads a = e/\^2TTe, which of course coincides with the result 
of Goldenfeld et all 



^It should be mentioned that thus obtained does not satisfy the Barlenblatt's equation even in 

the order of e; this is because the Gaussian form for F{x^ /t) is not correct. To obtain the precise form of 
F{x'^/t), one may insert Uj^{x,t) into the Barlenblatt's equation and solve F(^). 
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